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CALCULATION  OF  VELOCITY  PROFILES  AND  VISCOUS  HEATING  EFFECTS 
IN  PROPELLANT  EXTRUSION  DIES 

R.C.  Warren 


SUMMARY 


The  flow  of  a  propellant  dough  in  capillary  and  annular 
dies  has  been  analysed  with  a  finite  difference  ccagiuter 
program.  Velocity  profiles  have  been  calculated  and  the 
effects  of  viscous  heating  at  the  high  extrusion  rates 
used  in  propellant  processing  hav'  been  investigated. 
Viscous  heating  was  found  to  have  a  dramatic  effect, 
causing  considerable  rearrangeawnt  of  velocities  in  the 
die.  It  is  suggested  that  viscous  heating  reduces 
extrusion  pressures  and  die  swell  and  it  nay  have  a 
beneficial  effect  on  extrudate  surface  finish. 
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1.  INTRODUCTION 

Gun  propellants  can  be  laanufactured  by  a  range  of  processes,  smst  of  trtiicb 
involve  extrusion  throuf^  dies  to  produce  the  required  cross-sectional  shapes. 
Irregularities  in  diaensions  of  propellant  grains  caused  by  poor  extrusion  can 
result  in  unacceptable  variations  in  snizzle  velocity  and  gun  pressure.  Soae 
of  the  problens  encountered  in  extrusion  are  die  swell,  surface  roughness  and 
melt  fracture. 

Little  is  known  either  about  the  way  propellant  doughs  flow  in  dies,  or  about 
the  rheological  properties  of  propellant  doughs.  However,  they  SHist  have  a 
strong  Influence  on  extrusion  behaviour.  An  analysis  of  the  factors  affecting 
propellant  extrusion  is  necessary,  but  the  task  presents  considerable 
difficulties.  The  geometries  of  dies  are  often  co^licated  and  no 
satisfactory  constitutive  equations  describing  polymer  flow  are  yet  available. 
A  further  coBq>licating  factor  is  that  extrusion  often  takes  place  at  rates 
where  viscous  heat  generation  can  be  expected  to  occur.  In  the  extrusion  of 
polymers,  viscous  heating  is  known  to  cause  large  radial  temperature  gradients 
and  velocity  rearrangements (ref. 1).  A  comprehensive  review  of  viscous  heating 
has  been  given  by  Winter(ref .2) .  Quantitative  analysis  of  heating  effects  in 
dies  is  difficult,  as  viscous  heating  is  described  by  coupled,  non-linear 
differential  equations  which  require  numerical  awthods  for  solution. 

In  this  report,  a  finite  difference  program  has  been  used  to  calculate 
velocity  profiles  in  the  shape  forming  sections  of  three  propellant  extrusion 
dies.  The  dies  analysed  mre  a  solid  cord  die  for  casting  powder,  a  single  pin 
die  for  small  arms  propellant,  and  a  7  pin  die  for  cannon  propellant.  The 
effect  of  viscous  heating  in  the  solid  cord  and  single  pin  dies  was 
calculated.  Various  approxisMtions  were  made,  but  the  results  illustrate  the 
kinds  of  flow  which  may  be  encountered  in  extrusion,  and  the  importance  of 
viscous  heating. 


2.  PROGRAM 

The  mathematical  model  used  in  this  work  has  been  described  by  Cox  and 
Nacosko(ref . 1)  and  the  basic  computer  program  was  obtained  from  the  Ph.D 
thesis  of  Cox(ref.3).  The  model  includes  the  following  assumptions. 

(1)  The  flow  is  laminar  and  steady. 

(2)  The  flow  is  axisymmetric . 

(3)  There  is  no  fluid  slip  at  the  wall. 

(4)  All  normal  stresses  are  zero. 

(5)  Thermal  conductivity,  k,  heat  capacity,  Cp,  and  density,  p,  are 
independent  of  temperature  and  pressure. 

(6)  Viscosity,  q,  can  be  expressed  as  a  function  of  shear  rate,  y,  and 
tesqierature,  T. 


With  these  assumptions  the  equations  describing  flow  in  a  capillary  reduce  to 


8Vz  l8(r.Vr) 

—  4.  -  . —  s  0 

8z  r8r 


(1) 
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dP  l8(r.Trz) 

-  + - -  0 

(2) 

dz  rdr 

rVrai 

P^Cp  —  + 

-  8r 


IrzSVz  k8  (rdT) 

dr  r8z  8r 


(3) 


with  boundary  conditions:- 


—  =  trz  =  0  r  =  0 

Vz  =  Vr  =  0  r  =  R 

T  =  To  z  =  0 

Vz  =  VZfl  V  =  0 

Vr  =  0  z  =  0 

P  =  0  z  =  0 

8T  ^  -Hu(To-T)  r  =  R 

8r  R 


0  <  z  <  L 

0  <  z  <  0 
0  <  r  <  R 
0  <  r  <  R 
0  <  r  <  R 
0  <  r  <  R 
0  <  z  <  L 


R  and  L  are  the  radius  and  length  of  the  capillary,  To  the  initial  dough 
temperature,  Vz(r)  the  inlet  velocity  profile  assuming  isothermal  conditions 
at  the  inlet,  Nu  the  Nusselt  number  (2hR/k),  and  b  the  beat  transfer 
coefficient  at  the  wall.  A  complete  Notation  of  syid>ols  is  given  on  page  7. 
Cox  and  Hacosko  found  that  the  empirical  relationship 


Nu 


1.75(wCp/ia.) 


l/s 


1.75(Gz) 


i/s 


(4) 


could  be  used  to  specify  the  thermal  boundary  condition  at  the  wall.  6z  is 
the  Graetz  number  (wCp/kL)  and  w  the  awss  flow  rate(ref.l).  This  relationship 
leads  to  the  following  expression  for  heat  transfer  to  the  walls, 


k  i/s  (5) 

h  =  0.875— (  wCp/kL) 

R 


which  sllows  the  temperature  of  the  walls  to  increase  as  the  dough  tanperature 
increases . 

Equations  (1)  and  (2)  are  integrated  directly  by  nuMrical  methods,  and 
Equation  (3)  solved  by  an  iterative  procedure  after  converting  it  to  a  finite 
difference  equation.  The  differences  age  generated  on  a  variable  apaeiag  grid 
which  is  finest  where  rates  of  change  of  tiM  variables  are  higj^t.  An 


-  3  - 


l«itL-0283-TK 


initial  tenperature  distribution  is  assuaed  and  the  velocities  and  pressures 
calculated  from  equations  (1)  and  (2).  Equation  (3)  is  then  used  to  calculate 
a  new  tenperature  distribution  which  serves  as  the  basis  for  a  further 
iteration.  The  process  is  continued  until  satisfactory  convergence  is 
obtained. 

Cox's  original  prograu  was  written  for  capillary  flow  and  it  had  to  be 
extensively  nodified  to  enable  it  to  be  used  for  annular  flow  geoaetries.  In 
capillary  flow  the  integration  of  equation  (2)  is  straightforward.  In  annular 
flow  the  integration  has  to  be  split  into  two  parts,  one  part  over  the  radius 
between  the  inner  wall  and  the  radial  distance  of  the  uaxiaua  velocity, 
R(Vaux),  and  the  other  between  R(Vaax)  and  the  outer  wall.  RCVisax)  is  not 
known  a  priori  and  aiust  be  calculated  by  an  iterative  procedure  for  each  axial 
section. 


3.  DIE  GECMETRIES 

The  flow  path  of  dough  in  propellant  production  dies  is  quite  couplicated. 
Screens  and  sieves  are  used  to  collect  any  luiqis  or  ungelatinised  NC 
particles.  The  dough  then  flows  into  a  distribution  region  in  the  die  block 
containing  several  channels  leading  to  as  nany  as  64  separate  dies.  The  dies 
usually  taper  into  a  final  parallel  section  which  detecaiines  the  shape  of  the 
extrudate.  Dies  for  perforated  propellant  grains  have  pins  idiich  require 
coaq>licated  structures  to  support  then  and  these  cause  even  nore  coaiplex  flow 
patterns. 

The  flow  in  nonparallel  sections  of  a  die  is  intricate  and  cannot  be  readilly 
Bodelled  with  the  Cox  program.  However,  insight  into  flow  behaviour  and 
extrusion  problems  can  be  gained  by  studying  the  flow  in  the  die  parallel 
sections.  In  this  work  the  flow  at  the  entrance  of  the  die  parallel  section 
was  assumed  to  be  isothermal.  As  the  dough  progressed  through  the  die,  the 
velocity  and  tesperature  profiles  were  allowed  to  change  in  response  to 
viscous  heating. 

The  shapes  and  dimensions  of  the  dies  studied  are  illustrated  in  figure  1. 
The  flow  in  the  parallel  sections  of  dies  used  to  make  solid  cord  casting 
powder  for  CDB  propellants  and  in  dies  with  single  pins  used  to  make 
perforated  ssuill  arms  propellant  grains  can  be  readily  modelled.  However,  in 
multi-pin  dies,  such  as  those  used  for  multiperforsted  cannon  propellant,  the 
flow  is  not  axisymmetric  and  there  are  radial  and  circumferential  flows  caused 
by  the  outer  circle  of  pins.  To  obtain  a  picture  of  the  flow  in  this  type  of 
die,  flows  in  two  longitudinal  sections  of  the  die  were  studied.  One  section 
contained  a  diameter  intersecting  two  outer  pins,  section  A-A  in  figure  Ic, 
and  the  other  contained  a  diaawter  equidistant  from  the  outer  pins,  section  B- 
B  in  figure  Ic.  By  symmetry,  the  flows  in  these  two  sections  are  axial  and  can 
be  modelled  by  considering  them  to  be  sections  of  an  axisymmetric  die. 


4.  RHEOLOGICAL  DATA 

Measurements  of  the  rheological  propertiea  of  double  base  doughs  of  similar 
composition  to  casting  powder  doughs  have  been  made  previously(ref .4).  Shear 
stress,  surface  temperature  rise  and  die  swell_were  measured  at  shear  rates  at 
the  die  wall  ranging  from  from  1  s~^  to  1200  s~^.  The  Cox  program  was  tested 
by  using  the  shear  stress  data  to  predict  the  expected  surface  temperature 
rise  in  capillary  dies  of  various  lengths.  Good  agreement  was  obtained 
between  the  predicted  and  experimentally  observed  values,  as  is  illustrated  la 
figure  2. 

The  rheological  data  for  the  double  base  dough  wss  suitable  for  medelliag  flew 
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in  the  castint  powder  die.  However,  comparable  rheological  data  for  single 
and  triple  base  doughs  for  ■odelling  flows  in  aaall  aras  and  canaoa  powder 
dies  has  net  yet  been  deterained.  In  order  to  illustrate  the  type  of  flow 
«diich  occurs  in  these  dies,  the  data  for  the  double  base  dough  were  used. 


5.  RESULTS 


5.1  Casting  powder  die 

Calculated  velocity  profiles  at  the  die  inlet  and  exit,  and  the  teaperature 
profile  at  the  exit,  are  shown  in  figure  3.  The  dough  inlet  teaperature 
was  20®C.  Calculations  were  aade  for  an  extrudate  velocity  of  2.5  a/s,  the 
aaxiaua  value  used  in  factory  production.  The  effect  of  viscous  heating  on 
the  exit  profile  is  draawtic.  The  flow  is  plug-like,  with  virtually  a 
constant  velocity  for  95%  of  the  radius  of  the  die.  Alaost  all  of  the 
change  of  velocity  occurs  in  the  outer  5%  of  the  radius  near  the  die  wall. 
The  teaperature  rise  is  largely  confined  to  the  outer  10%  of  the  radius, 
and  reaches  50^C  at  the  wall.  The  flow  rate  is  so  high  there  is 
insufficient  tiae  for  significant  heat  conduction  to  the  interior  of  the 
dough. 

5.2  Saall  aras  propellant  die 

A  flow  rate  corresponding  to  a  typical  production  extrudate  velocity  of 
0.8  a/s  was  assuaed.  The  velocity  profiles  at  the  inlet  and  exit,  and  the 
teaperature  profile  at  the  exit,  are  shown  in  figure  4  for  a  dough  inlet 
teaperature  of  20^C.  The  shear  rate  is  given  by  the  slope  of  the  velocity 
profile  and  it  can  be  seen  that  the  shear  rate  is  auch  greater  at  the  pin 
than  at  the  outer  wall  at  both  the  die  inlet  and  exit.  The  effect  of 
viscous  heating  is  to  reduce  the  shear  rate  at  the  centre  of  flow  and  to 
increase  the  shear  rate  at  the  walls,  particularly  at  the  pin.  The 
teaperature  rises  are  substantial,  46**C  at  the  pin,  and  30^C  at  the  outer 
wall. 

5.3  Cannon  powder  die 

The  flow  rates  used  for  the  coaaercial  extrusion  of  cannon  propellant  are 
considerably  lower  than  those  used  for  saall  anas  propellants  and  for 
casting  powder.  Consequently  the  effect  of  viscous  heating  can  be  ignored. 
A  flow  rate  corresponding  to  an  extrudated  velocity  of  0.025  a/s  was  used 
in  the  analysis. 

The  velocity  profiles,  assuaing  an  inlet  teaperature  of  20^  in  the  two 
sections  of  the  die  shown  in  figure  1,  are  illustrated  in  figure  5.  The 
shear  rate  is  low  near  the  outer  wall,  but  near  the  pins  it  is  considerably 
higher.  The  velocity  profiles  near  the  inner  and  outer  surfaces  of  the 
outer  ring  of  pins  are  different  which  indicates  that  flow  near  the  outer 
pins  is  not  synnetric.  This  would  tend  to  cause  the  outer  holes  in  the 
extrudate  to  be  non-circluar. 


6.  DISCUSSIOH 

Although  the  calculated  velocity  and  teaperature  profiles  can  only  illustrate 
the  patterns  of  flow  expected  to  occur  in  typical  propellant  dies,  seas 
conclusions  can  be  drawn. 

Viscous  heating  has  a  very  strong  effect  mi  propellant  flow.  The  increase  in 
teaperature  causes  a  reduction  in  extrusion  stresses,  vhic)i  in  tun  reduesa 
the  pressure  required  to  extrude  at  a  given  rate.  This  is  illastntad  in 
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figure  2,  taken  froa  reference  4,  which  ahows  observed  and  calculated  shear 
stresses  vs  shear  rate  for  a  double  base  dough  extruded  through  a  capillary 
die. 

The  low  viscosity  of  the  surface  layers  of  extrudates  due  to  viscous  heating 
should  contribute  to  an  iag»roveawnt  in  the  surface  finish  of  propellant 
grains.  Surface  finish  of  polyner  extrudates  is  often  swrred  by  'sharkskin', 
fine  regular  fissures  in  the  surface  which  occur  %dien  the  extrusion  rate 
exceeds  a  critical  value.  Cogswell  has  shown  that  sharkskin  is  due  to 
stretching  of  the  surface  layer  of  the  extrudate  as  it  leaves  the  die(ref.5). 
At  the  lip  of  the  die  the  velocity  of  the  surface  of  the  extrudate  goes  froa 
near  zero  to  a  relatively  high  value  in  a  snail  distance.  This  acceleration 
of  flow  causes  large  extensional  stresses  in  the  outer  layers  of  the 
extrudate.  Vdien  the  extensional  stresses  exceed  the  strength  of  the  material 
fracture  of  the  surface  occurs.  However,  if  the  surface  layer  has  a  low 
viscosity  the  extensional  stresses  will  be  reduced  and  the  likelihood  of 
sharkskin  occur ing  is  consequently  aininised. 

The  variation  of  die  swell  with  shear  rate  for  the  double  base  dough  extruded 
though  a  capillary  die  is  illustrated  in  figure  6(ref.4).  For  each  die  ,  the 
swell  increases  with  rate  at  low  shear  rates,  but  it  peaks  and  then  decreases 
at  high  shear  rates.  Longer  dies  show  a  greater  decrease  in  swell  at  high 
rates.  The  shear  rates  at  which  die  swell  peaks  correspond  to  the  rates  at 
which  viscous  heating  becomes  important.  This  suggests  that  the  effect  of 
viscous  heating  is  to  reduce  die  swell. 

There  are  two  possible  ways  that  viscous  heating  can  affect  die  swell.  Firstly 
by  changing  the  viscosity  gradient  across  the  die,  and  secondly  by  changing 
the  relaxation  tines  of  the  elastic  stresses. 

Phuoc  and  Tanner  used  a  finite  element  program  to  calculate  die  swell  for  a 
Newtonian  fluid  with  a  temperature  dependent  viscosity  flowing  in  a  die  with 
constant  wall  tenperature(ref .6) .  Under  these  circumstances  the  centre  of  the 
fluid  becomes  hotter  than  the  walls  and  there  is  an  increasing  viscosity 
gradient  from  the  centre  to  the  walls.  Their  calculations  showed  that  die 
swell  would  increase  with  increasing  shear  rate  as  a  consequence  of  the 
viscosity  gradient  across  the  die.  For  the  cases  considered  in  this  study, 
the  temperature  gradient  is  opposite  to  that  of  Phuoc  and  Tanner,  that  is,  the 
viscosity  of  the  outer  layer  is  less  than  in  the  interior.  Under  these 
circumstances  the  effect  of  viscous  beating  would  be  to  reduce  die  swell. 

Die  swell  has  been  related  to  the  relaxation  of  elastic  stresses  induced  in 
the  extrudate  by  shear  flow  in  the  die(ref.7).  Viscous  heating  tends  to  cause 
plug  flow  where  most  of  the  extrudate  travels  as  an  almost  rigid  body  which 
would  not  give  rise  to  die  swell.  Large  shear  stresses  are  confined  to  the 
ssmII  volume  of  low  viscosity  amterial  near  the  wall.  They  should  be  able  to 
relax  quickly  because  of  the  liigh  temperature  and  so  reduce  the  amount  of  die 
swell . 

The  fact  that  teaq>erature  rises  are  confined  to  areas  near  the  die  walls  and 
pins  indicates  that  the  surface  of  the  extrudate  may  have  properties  different 
from  the  bulk  material.  The  very  high  shear  rates  at  the  surface  would  tend 
to  cause  molecular  orientation  at  the  surface  which  would  be  frozen  in  as  the 
extrudate  cooled.  It  is  also  likely  that  the  degree  of  gelatinisation  of  a 
propellant  is  affected  by  extrusion  tes^erature,  which  should  result  in  the 
surface  of  the  extrudate  being  more  highly  gelatinised  then  the  interior. 

Sections  of  solid  cord  extruded  at  300  s~^  when  viewed  under  a  microscope  with 
crossed  polaroids  show  a  bright  layer  at  the  surface,  see  figure  7.  This 
supports  the  assumption  that  molecular  orientation  occurs  at  the  surfiM 
during  extrusion  of  propellant  doughs. 
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The  oriented  layer  «t  the  surface  of  the  propellant  grains  could  be  expected 
to  affect  the  diffusion  of  deterrents »  and  also  ispact  fracture  toughness. 
Further  research  is  necessary  to  ascertain  whether  or  not  these  effects  are 
significant. 

The  largest  source  of  error  in  the  calculations  in  this  report  is  the 
assuaption  of  the  rheological  properties  of  the  dough.  It  is  necessary  to 
detensine  the  viscosity  of  the  propellants  doughs  of  interest  under  the 
conditions  which  occur  in  production  for  accurate  calculations  to  be  aade. 
However,  the  results  obtained  suggest  that  it  aay  be  possible  to  design  dies 
and  to  vary  the  rheological  properties  of  the  dough  so  that  viscous  heating 
could  be  used  to  iaprove  propellant  quality  by  reducing  die  swell  and  surface 
roughness . 


7.  CONCLUSICmS 

Velocity  and  teaperature  profiles  of  propellant  flow  in  three  typical 
propellant  production  dies  have  been  calculated.  It  has  been  found  that 
viscous  heating  causes  large  teaperature  rises  at  the  walls  of  the  die  and 
near  pins  and  that  these  rises  have  a  significant  effect  on  flow  behaviour. 

It  is  suggested  that  viscous  heating  could  be  used  to  iaprove  surface  finish 
and  to  reduce  die  swell  of  extrudates.  However,  considerably  aore  experiaental 
and  theoretical  work  is  required  before  the  controlled  use  of  viscous  heating 
could  be  utilised  as  an  aid  in  propellant  processing. 
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NOTATION 

heat  capacity 

Graetz  rniaber 

heat  transfer  coefficient 

therml  conductivity 

die  length 

Nusselt  nuaber 

pressure 

radial  coordinate 
die  radius 
temperature 
inlet  teaqperature 
axial  velocity 
radial  velocity 
mass  flow  rate 
axial  coordinate 
shear  rate 
viscosity 
density 
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(a)  Casting  powler  die 


(b)  Snallanns  propellant  die 


(c)  Cannon  powder 


Figure  1.  Diaensions  and  geoaetries  of  dies 
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Apparent  shear  rate  (s  ”* ) 

Figure  2.  Shear  stress  as  apparent  shear  rate(ref.4) 
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